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Abstract 

We analyze the chargino contributions to Bd — Bd mixing and CP asymmetry of the 
Bd — » J/ipKs decay, in the framework of the mass insertion approximation. We derive 
model independent bounds on the relevant mass insertions. Moreover, we study these 
contributions in supersymmetric models with minimal flavor violation, Hermitian flavor 
O ■ structure, and small CP violating phases and universal strength Yukawa couplings. We 

show that in supersymmetric models with large flavor mixing, the observed values of 
sin 2(3 may be entirely due to the chargino-up-squark loops. 



1 Introduction 

X 

Since its discovery in 1964 in the A"-meson decays, the origin of CP violation remains an 
open question in particle physics. In standard model (SM), the phase of the Cabbibo- 
Kobayashi-Maskawa (CKM) quark mixing matrix provides an explanation of the CP 
violating effect in these decays. Although the SM is able to account for the observed CP 
violation in the kaon system and the recent measurement of the (time-dependent) CP 
asymmetry in Bd — > J/ipKs decays, new CP violating sources are necessarily required 
to describe the observed baryon asymmetry [1]. Moreover, it is expected that with B 
factories, the .B-system will represent an ideal framework for crucial tests of the CP 
violation in SM and probing new physics effects at low energy. 

It is a common feature for any new physics beyond the SM to possess additional CP 
violating phases beside the Sqkm phase. In supersymmetric (SUSY) models, the soft 
SUSY breaking terms contain several parameters that may be complex, as may also the 
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SUSY preserving \i parameter. These new phases have significant implications for the 
electric dipole moment (EDM) of electron, neutron, and mercury atom [2]. It was shown 
that the EDM can be suppressed in SUSY models with small CP phases [3, 4] or in SUSY 
models with flavor off-diagonal CP violation [3,5]. 

The idea of having small CP phases ( < 10~ 2 ) as an approximate CP-symmetry at low 
energy, could be an interesting possibility if supported by a mechanism of CP-symmetry 
restoration at high energy scale. However, this mechanism might also imply that the 
&ckm phase is small [6]. The large asymmetry of the £?-meson decay aj/^K s observed by 
BaBar and Belle [7] are in agreement with SM predictions for large Sckm, and thus the 
idea of small phases might be disfavored. 

However, in Ref. [8] it was shown that in the framework of the minimal supersymmetric 
standard model (MSSM) with non-universal soft terms and large flavor mixing in the 
Yukawa, supersymmetry can give the leading contribution to aj/^K s with simultaneous 
account for the experimental results in i^-system. Thus the supersymmetric models with 
small CP violating phases at high energy scale are still phenomenological viable. An 
alternative possibility for suppressing the EDMs is that SUSY CP phases has a flavor 
off-diagonal character as in the SM [5,9]. Such models would allow for phases of order 
(9(1) which may have significant effect in 5-physics [5]. 

A useful tool for analyzing SUSY contributions to flavour changing neutral current 
processes (FCNC) is provided, as known, by the mass insertion method [10]. One chooses 
a basis for the fermion and sfermion states where all the couplings of these particles to 
neutral gauginos are flavour diagonal, leaving all the sources of FC inside the off-diagonal 
terms of sfermion mass matrix. These terms are denoted by (A AB y j , where as usual 
A, B = (L, R) and i,j — 1, 3 indicate chiral and flavour indices respectively and q — u,d. 
The sfermion propagator is then expanded as a series of (5 AB )ij = (A AB y j /rh 2 , where m 2 
is an average sfermion mass. This method allows to parametrize, in a model independent 
way, the main sources of flavour violations in SUSY models. In this framework, the gluino 
and chargino contributions to the i^-system have been analyzed in references [10] and 
[11] respectively. These analyses showed that the bounds on imaginary parts of mass 
insertions, coming from gluino exchanges to Ek and e'/e, are very severe [10], while the 
corresponding ones from chargino exchanges are less constrained [11]. In particular, in or- 
der to saturate Ek from the gluino contributions one should have [10] \J ' \lm(5i 2 ) LL \ ~ 10~ 3 

or \J \Im(5 1 2 ) lr\ ~ 10 -4 , and y / |Im(5f 2 )i J R| ~ 10~ 5 from e'/e, while chargino contributions 
require [11] Im(<5" 2 )| L ~ 10~ 2 , for average squark masses of order of 500 GeV and gluino 
masses of the same order. 

Recently, in the framework of mass insertion approximation, gluino contributions to 
the Bd — Bd mixing and CP asymmetry in the decay Ba — > J / ipKs have been analyzed by 
including next-to-leading order (NLO) QCD corrections [12] (see also Ref. [13]). However, 
an analogous study for chargino contributions to these processes is still missing. This kind 
of analysis would be interesting for the following reasons. First, it would provide a new 
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set of upper bounds on the mass insertion parameters, namely (S^)ab in the up-squark 
sector, which are complementary to the ones obtained from gluino exchanges (which only 
constrain (<^-)ab)- Second, upper bounds on (<%)ab would be very useful in order to 
perform easy tests on SUSY models which receive from chargino exchanges the main 
contributions to Bd — Bd and CP asymmetry Indeed, in many SUSY scenarios the gluino 
exchanges are always sub-leading. 

In this paper we focus on the dominant chargino contributions to the Bd — Bd mixing 
and CP asymmetry aj/^K s - We use the mass insertion method and derive the correspond- 
ing bounds on the relevant mass insertion parameters. We perform this analysis at the 
NLO accuracy in QCD by using the results available in Ref.[12]. As an application of our 
analysis, we also provide a comparative study for supersymmetric models with minimal 
flavor violation, Hermitian flavor structure with small CP violating phases and univer- 
sal strength of Yukawa couplings. We show that in all these scenarios, by comparing 
(<%)ab and (5^)ab with their corresponding upper bounds, the chargino contributions 
are dominant over the gluino ones. 

The paper is organized as follows. In section 2 we present the supersymmetric con- 
tributions to B d — Bd mixing and CP asymmetry aj/^Ks- We start with a brief review 
on gluino contributions and then we present our results for the chargino ones, both in 
the mass insertion approach. In section 3 we derive model independent bounds on the 
relevant mass insertions involved in the Bd — Bd mixing and aj/^x s - In section 4 we 
generalize these results by including the case of a light stop-right. Section 5 is devoted to 
the study of the supersymmetric contribution to aj/^Ks * n three different supersymmetric 
models. We show that the observed values of sin 2f3 may be entirely due to the chargino— 
up-squark loops in some classes of these models. Our conclusions are presented in section 
6. 



2 Supersymmetric contributions to AB = 2 transitions 

We start this section by summarizing the main results on B d — B d mixing and CP asym- 
metry aj/.,pK s , then we will consider the relevant SUSY contributions to the effective 
Hamiltonian for AB = 2 transitions, given by the chargino and gluino box diagram 
exchanges. 

In the Bd and Bd system, the flavour eigenstates are given by Bd = (bd) and Bd = 
(bd). It is customary to denote the corresponding mass eigenstates by Bh = pBd + qBd 
and B L = pB d — qBd where ind ices H and L refer to heavy and light m ass eigenstates 
respectively, and p — (1 + e B )/ ^2(1 + \e B \), q = (1 — e B )/^2(l + \e B \) where e B is the 
corresponding CP violating parameter in the Bd — Bd system, analogous of e in the kaon 
system [14]. Then the strength of B d — B d mixing is described by the mass difference 



AM Bd = M Bh -M Bl . 
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whose present experimental value is AM^ = 0.484 ± 0.010 (ps) 1 [14]. 

The CP asymmetry of the B d and B d meson decay to the CP eigenstate 4>K S is given 

by 

= T(m)->1>Ks) + T(B d (t)^il>Ks) = * KS Sm(Am ^ t} - (2) 
The most recent measurements of this asymmetry are given by [7] 

Oj, Kg = 0.59 ±0.14 ±0.05 (BaBar) , 

a^ Ks = 0.99 ±0.14 ±0.06 (Belle). (3) 

where the second and third numbers correspond to statistic and systematic errors respec- 
tively, and so the present world average is given by a^K s = 0-79 ± 12. These results 
show that there is a large CP asymmetry in the B meson system. This implies that 
either the CP is not an approximate symmetry in nature and that the CKM mechanism 
is the dominant source of CP violation [5] or CP is an approximate symmetry with large 
flavour structure beyond the standard CKM matrix [8]. Generally, AM Bd and cl^k s can 
be calculated via 

AM Bd = 2\(B° d \H^\B° d )\ , (4) 
a^ Ks = sin2& ff , and fa = ^Tg(B° d \H^ B=2 \B d ) . (5) 

where H^ B=2 is the effective Hamiltonian responsible of the AB = 2 transitions. In the 
framework of the standard model (SM), cl^k s can be easily related to one of the inner 
angles of the unitarity triangles and parametrized by the Vckm elements as follows 



« s = sin2/3, /3 = ar g(-^|], (6) 



In supersymmetric theories the effective Hamiltonian for AB = 2 transitions, can 
be generated, in addition to the W box diagrams of SM, through other box diagrams 
mediated by charged Higgs, neutralino, photino, gluino, and chargino exchanges. The 
Higgs contributions are suppressed by the quark masses and can be neglected. The 
neutralino and photino exchange diagrams are also very suppressed compared to the gluino 
and chargino ones, due to their electroweak neutral couplings to fermion and sfermions. 
Thus, the dominant SUSY contributions to the off diagonal entry in the 5-meson mass 
matrix, M 12 (B d ) = (B° d \H^ B = 2 \B° d ), is given by 

M 12 (B d ) = Mff(B d ) + MUBd) + M\ 2 + {B d ). (7) 

where Mf^(B d ), M.\ 2 {B d ), and M.\ 2 indicate the SM, gluino, and chargino contributions 
respectively. The SM contribution is known at NLO accuracy in QCD [14] (as well as the 
leading SUSY contributions [12]) and it is given by 



Mff(B d ) = ^ B B B JlM Bd M 2 w {V td V; b fS,{x t ), (8) 
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where fs d is the B meson decay constant, B^ d is the renormalization group invariant B 
parameter (for its definition and numerical value, see Ref. [14] and reference therein) and 
r) = 0.55 ± 0.01. The function S (x t ), connected to the AB = 2 box diagram with W 
exchange, is given by 

Ax t - \\x\ + xf 3xf\nx t , 

* ( *> = 4(i -x.)' — (9) 

where x t = Mf/M^. 

The effect of supersymmetry can be simply described by a dimensionless parameter 
r\ and a phase 26 d defined as follows 

r 2 r 2i9 d _ M 12 (B d ) 

where AM Bd = 2\Mff(B d )\rj. Thus, in the presence of SUSY contributions, the CP 
asymmetry B d — > ipK s is modified, and now we have 

a fKs = sin2/3 cff = sin(2/3 + 29 d ) . (11) 

Therefore, the measurement of cl^k s would not determine sin 2f3 but rather sin 2/3 e g , where 



Mf™ Y (B d y 

Mff(B d ) 



29 d = arg ( 1 + a " sm /d? ) » ( 12 ) 



and M? 2 USY (£ d ) = A4f 2 (B d ) + Mt(B d ). 



2.1 Gluino contributions 

The most general effective Hamiltonian for AB = 2 processes, induced by gluino and 
chargino exchanges through AB = 2 box diagrams, can be expressed as 

#e A / =2 = E^WiM +^Ci(fi)Qi(fi) + h.c. , (13) 

i=i i=i 

where Cj(/i), Ci(fi) and Qi(fi), Qi(^) are the Wilson coefficients and operators respectively 
renormalized at the scale /i, with 

Q l = d a L ^bl dh.bl Q 2 = d a R b a L 4&£, Q 3 = d a R b{ d£&£, 

Q 4 = d a R b a L 4b p R , Q 5 = d a R bUlb%. (14) 

In addition, the operators <5i,2,3 are obtained from <5i,2,3 by exchanging L <-> /?. 

Now we summarize the main results for gluino contributions to the above Wilson 
coefficients at SUSY scale, in the framework of mass insertion approximation. As we 
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will show in the next section, in order to connect the Wilson coefficients at SUSY scale 
with the corresponding ones at low energy scale fx ~ 0(rrib), the RGE equations for QCD 
corrections must be solved. In the case of the gluino exchange all the above operators 
give significant contributions and the corresponding Wilson coefficients are given by [10], 
[12] 



Ci(M s ) = -- j |- ? (24x/ 6 (x) + 66/ 6 (x))(^ 3 )| L , 



216m? 
n 2 

° 3{Ms) = ^^ 36 ^)(<^)k, (15) 
C,(M S ) = -^-2 {(504xf 6 (x) - 72/ 6 (x)) (5 d 13 ) LL (5 d 13 ) RR - 132f 6 (x)(5 d l3 ) LR (5 d l3 ) RL 

C 5 (M S ) = [(2Axf e (x) + 120/ 6 (x)) (5 d 3 ) LL (5 d 3 ) RR - 180Ux)(5 d 3 ) LR (5 d 3 ) RL 

where x = m^/m^ and m 2 is an average squark mass. The expression for the functions 
fe(x) and fe(x) can be found in Ref.[12]. The Wilson coefficients C\s are simply obtained 
by interchanging L <-> R in the mass insertions appearing in Ci_3. 



2.2 Chargino contributions 

Here we present our results for the chargino contributions to the effective Hamiltonian 
in Eq.(13) in the mass insertion approximation. The leading diagrams are illustrated in 
Fig. 1, where the cross in the middle of the squark propagator represents a single mass 
insertion. As we will explain in more details below, the dominant chargino exchange, can 
significantly affect the operators Qi and Q3 only. We remind here that in the case of the 
K — K mixing, the relevant chargino exchange affects only the operator Qi [11], as in the 
SM. 

In the framework of mass insertion approximation, one chooses a basis (super-CKM 
basis) where the couplings of the fermions and sfermions to neutral gauginos are flavour 
diagonal. In this basis, the interacting Lagrangian involving charginos is given by 

Asx+ = -9 EE Oi K la dl (X + T u b L - U* k2 {Yt g .K + ) ab 4 (x + T u\ 

k a,b 

-V: 2 (K.Y u dia % b d a L (x + Yu b R ). (16) 

where Y^ g are the diagonal Yukawa matrices, and K is the usual CKM matrix. The 
indices a, b and k label flavour and chargino mass eigenstates respectively, and V, U are 
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Figure 1: The leading chargino-up-squark contribution to the Bj — Bj mixing. 



the chargino mixing matrices defined by 

U*M^+V~ l = diag(m~+,m-+), and = 



M 2 V2M W sin (5 

V2M W cos (3 (j, 



(17) 



As one can see from Eq.(16), the Higgsino couplings are suppressed by Yukawas of the light 
quarks, and therefore they are negligible, except for the stop-bottom interaction which is 
directly enhanced by the top Yukawa (Yt)- The other vertex involving the down and stop 
could also be enhanced by Y t , but one should pay the price of a A 3 suppression, where 
A is the Cabibbo mixing. Since in our analysis we adopt the approximation of retaining 
only terms proportional to order A, we will neglect the effect of this vertex. Moreover, 
we also set to zero the Higgsino contributions proportional to Yukawa couplings of light 
quarks with the exception of the bottom Yukawa Y b , since its effect could be enhanced 
by large tan/5. In this respect, it is clear that the chargino contribution to the Wilson 
coefficients C4 and C5 is negligible. Furthermore, due to the colour structure of chargino 
box diagrams there is no contribution to C2 or C2. However, as we will show in the next 
section, chargino contributions to C2 or C2 will be always induced at low energy by QCD 
corrections through the mixing with C 3 . 

Now we calculate the relevant Wilson coefficients Ci 3 (M s ) at SUSY scale M s , by 
using the mass insertion approximation. As mentioned in this case the flavor mixing 
is displayed by the non-diagonal entries of the sfermion mass matrices. Denoting by 
A q AB the off-diagonal terms in the sfermion (q = u, d) mass matrices for the up, down 
respectively, where A,B indicate chirality couplings to fermions A,B = (L,R), the A-B 
squark propagator can be expanded as 



<ft?t> = i (k 2 l - rhH - A^)" 1 



k 2 — fh? (k 2 



ABJab 



m 2 ) 2 



+ 0(A 2 ), (18) 



where q = u,d selects up, down sector respectively, a,b — (1, 2, 3) are flavour indices, 1 
is the unit matrix, and fh is the average squark mass. As we will see in the following, 
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it is convenient to parametrize this expansion in terms of the dimensionless quantity 
{5 q AB )ab = (A AB ) >ab I 'm 2 ■ At the first order in the mass insertion approximation, we find 
for the Wilson coefficients Ci 3 (Ms) the following result 

CUMs) = ^^E{l^i| 2 |^i| 2 ((^ i )i 1 + 2A(^)3i(^L)3 2 ) 



m^V^V* i(5 LL ) 31 (5 u RL ) 31 + A(^ L ) 32 (^ L ) 31 + A(^) 3 i(^ L : 



32 



+ Y'VaV^V*^^ + 2A(^ L ) 31 (<^ L ) 32 ) }L 2 (x i ,x j ), (19) 



C ^ MS) = Tq^^^^ Ui2Uj2VjlVil i^LL)ll + 2 H^LL)3l( S LL)32)Lo{x i ,X j ), (20) 



1927T 2 m 2 



hJ 



where Xi = m^ + /m 2 , and the functions L (x,y) and L 2 (x,y) are given by 

— f x h (x) - y h (y)\ . N -11 + 7a; - 2a; 2 61nx 
M*,v) = — J, Mx)= (1 _ x)3 

x h 2 (x) - y h 2 (y) , . 2 + 5a; - a; 2 6a; In a; , . 

L2(X ' y) = a^ > W= (l-af + (T^ ^ 

As in the gluino case, the corresponding results for C\ and C\ coefficients are simply 
obtained by interchanging L <-> i? in the mass insertions appearing in the expressions for 

Cl,3- 



3 Constraints from AMs d and sin 2/3 

In this section we present our numerical results for the bounds on (S AB )ij which come 
from AM Bd and CP violating parameter sin 2f3. We start with the chargino contributions 
which is found to be the dominant SUSY source in various models [5,8,15]. We also 
provide analytical expressions for AM Bd and sin 2f3 as functions of the mass insertions in 
the Wilson coefficients Ci(M s ) of Eqs. (19), (20). In our calculation we take into account 
NLO QCD corrections in both Wilson coefficients and hadronic matrix elements given in 
[12]. 

In order to connect Cj(Ms) at SUSY scale Ms with the corresponding low energy ones 
Ci(n) (where \x ~ Oim^)\ one has to solve the renormalization group equations (RGE) 
for the Wilson coefficients corresponding to the effective Hamiltonian in (13). Then, Cj(/i) 
will be related to C;(M S ) by [12] 

CM = E E + VcP s) ) V a *C s (M s ), (22) 

i s 

where Ms > m t and r] = as{Ms)/as{n). The values of the coefficients bf ,s \ cf' s \ and 
appearing in (22) can be found in Ref. [12]. In our analysis the SUSY scale, where SUSY 
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particles are simultaneously integrated out, is identified with the average squark mass rh. 
By using the NLO results of [12], we obtain, for the relevant chargino contributions 

dip) = xM C^Ms), CM = x 2 {n) C 3 (M S ), CM = xM C 3 (M S ), (23) 

while for the other coefficients Cj(/i) = 0, (i — 4, 5). Numerical values for Xi(fx), evaluated 
at /j, — mb, are shown in table (1) for some representative values of Ms- Notice that 



the coefficients bf^ and cf a> are different from zero and so the contribution to C 2 (/jl) is 
radiatively generated at NLO by the off-diagonal mixing with C 3 (Ms). For the coefficients 
Ci_ 3 hold the same results as in Eq.(23) and in table (1), since the corresponding b\ r ' s ^ 
and cf"'^ coefficients in Eq.(22) are the same as the ones for the evolutions of Ci_ 3 [12]. 
The off diagonal matrix elements of the operators Qi are given by [12] 



.(23) 



{Bd\Qi\B d ) 
(B d \Q 2 \B d ) 

(B d \Q 3 \B d ) 

(B d \Q 4 \B d ) 

(B d \Q 5 \B d ) 



: m B d f Bd B M, 



24 \m h (n) + m d (/jL) / 

j_ ( m B d y 

24 \m h (n) + m d (/i) / 



m B J B B 2 (n), 



m B d f B BM, 



(24) 




12 \m b (n) + m d (fi) 



m B d f Bd BM). 



The value of B\ has been extensively studied on the lattice [16], but for the other Bi 
parameters, they have been recently calculated on the lattice by the collaboration in Ref. 
[17]. In our analysis we will use the following central values reported in [12], namely 
Bi(n) = 0.87, B 2 (n) = 0.82, B 3 {fi) = 1.02, B A (n) = 1.16, and B 5 = 1.91. The same 



M s 


Xl(fJ>) 


x 2 (/i) 


x 3 (/i) 


200 


0.844 


-0.327 


0.571 


400 


0.827 


-0.367 


0.536 


600 


0.817 


-0.389 


0.518 


800 


0.810 


-0.404 


0.506 



Table 1: Numerical values for the coefficients Xi (with i = 1,2,3) in Eq.(23) for some 
representative values of SUSY scale M s , and evaluated at the low energy scale jj, = m b 

results of Eq.(24) are also valid for the corresponding operators Qi, with same values for 
the Bi parameters, since strong interactions preserve parity. 
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Now we start our analysis, by discussing first the dominant chargino contribution to 
AM Bd . Using Eqs.(4), (19)-(20), and (23)-(24), we obtain for AM Bd the following result 

R = ((^l)3i + 2A(^ l )3i(^) 32 )x 

(2A 1 x 1 (At) J B 1 ( A t) + A 4 X(pi) (x 3 (ji)B 3 (ji) - hx 2 {fi)B 2 {fji)) ) 

+ ((«L)S1^)31 + A((^ £ )si(^)s2 + (^)32(^ L ) 3 l)) 2A 2 X 1 (|i)S 1 (|i) 

+ ((^l)3i + 2A(^ l )3i(^ l ) 32 ) 2A 3 x 1 (/i)B 1 (/i) (26) 



where R, which parametrizes the contributions of Qi- 3 operators, is obtained from R 
by exchanging L <-> _R in the mass insertions, = (m Bd /(m b (fi) + m d (fi))) 2 , and the 

expressions of are given by 

A 3 = Y^VaV^Vj^L^Xj), A 4 = Y?'£U i < l U j2 V jl V il L (x i ,x j ) (27) 



where the definition of the quantities appearing in (27) can be found in section [2]. Notice 
that the renormalization scheme dependence in Eq.(25) (for \x varying in the range fi ~ 
(m b /2,2m b )), is strongly reduced due to the NLO QCD accuracy. 

As customary in this kind of analysis [10], in order to find conservative upper bounds 
on mass insertions, the SM contribution to AM Bd is set to zero. Moreover, since we are 
analyzing AMg d which is a CP conserving quantity, we keep the squark mass matrices 
real. Upper bounds are then obtained by requiring that the contribution of the real 
part of each independent combination of mass insertions in Eq.(25) does not exceed the 
experimental central value AM Bd < 0.484(ps) _1 . 

These constraints depend on the relevant MSSM low energy parameters, in particular, 
by m, M 2 , fx and tan/9.*. Notice that with respect to the gluino mediated FCNC pro- 
cesses, which are parametrized by m, M 3 , the chargino mediated ones contains two free 
parameters more. 

In tables (2) and (3), we present our results for upper bounds on mass insertions 
coming from AM Bd , given for some representative values of m and M 2 and for fixed values 



of fi = 200 GeV and tan (3 = 5. In table (2) we provide constraints on y |Re [(5^) 3 i] | for 
several combinations of m and M 2 . We find that these bounds are almost insensitive to 
fx and tan (3 in the ranges of 200 — 500 GeV and 3-40 respectively. This can be simply 
understood by noticing that the contributions to AB = 2 transitions mediated by LL 



*With abuse of notation, we used here the same symbol fi for the renormalization scale of Wilson 
coefficients and the Higgs mixing parameter of MSSM. 
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M 2 \ m 


300 


500 


700 


900 


150 


1.3X10- 1 


1.7X10- 1 


2.2X10- 1 


2.8X10- 1 


250 


1.9X10" 1 


2.3X10" 1 


2.7X10" 1 


3.2X10" 1 


350 


2.7X10- 1 


2.8X10- 1 


3.3X10- 1 


3.7X10- 1 


450 


3.6X10" 1 


3.6X10" 1 


3.9X10" 1 


4.3X10" 1 



Table 2: Upper bounds on y^Re [(^J 3 i] 2 | from AM Bd (assuming zero CKM and SUSY 
phases), for \i = 200 GeV and tan (3 = 5, and for some values of rh and M 2 (in GeV). 

interactions are mainly given by the weak gaugino component of chargino. Therefore, the 
corresponding bounds are more sensitive to M 2 instead of \x and tan/3, since these last 
two parameters contribute to the Higgsino components of chargino. The only term in 
Eq.(25) which is quite sensitive to tan/5 is At, because it is proportional to the bottom 
Yukawa coupling squared. However, (S^ L ) 3 i, in addition to A±, receives contributions also 
from the A 1 term. This term is larger than A 4 and almost insensitive to tan/3, leaving 
the bounds on (6^)31 almost independent from tan/3. 



m 


\/Re[fe)3i] 2 


\l Re W RL hi? 


yj Re [(^ L )3l(^)32] | 


200 


1.4X10- 1 


4.7X10- 1 


2.1 xlO- 1 


400 


1.8X10- 1 


9.0X10- 1 


2.7X10" 1 


600 


2.2X10- 1 


1.5 


3.4X10" 1 


800 


2.7X10" 1 


2.3 


4.1X10" 1 



Table 3: Upper bounds on mass insertions as in table (2), for M 2 — /i — 200 GeV and 
tan (3 — 5. 



m 


>/|ne[(^L)si(^)si]| 


Vl^ [(^31(^)32] | 




200 


1.8X10- 1 


4.0X10" 1 


7.1 xlO" 1 


400 


3.0X10- 1 


6.3X10- 1 


1.3 


600 


4.5X10" 1 


9.5X10- 1 


2.3 


800 


6.3X10" 1 


1.3 


3.5 



Table 4: Upper bounds on mass insertions as in table (2), for M 2 = /i = 200 GeV and 
tan/3 = 5. 

In tables (3) and (4) we give our results for the real parts of the other mass insertions 
(and as well as for J\Re [(Sll)3i] 2 |) which are less constrained, for several values of rh 
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M 2 \ m 


300 


500 


700 


900 


150 


1.5X10- 1 


2.0X10- 1 


2.6X10- 1 


3.1X10" 1 


250 


2.2X10" 1 


2.6X10" 1 


3.1X10" 1 


3.6X10- 1 


350 


3.0X10- 1 


3.3X10- 1 


3.7X10- 1 


4.2X10- 1 


450 


4.0X10" 1 


4.1X10" 1 


4.4X10" 1 


4.8X10" 1 



Table 5: Upper bounds on y |Im [(5^ L ) 3 i] 2 | from sin2/3 = 0.79 (assuming a zero CKM 
phase), for /i = 200 GeV and tan/3 = 5, and for some values of m and M 2 (in GeV). 

and evaluated at M 2 — /i — 200 GeV and tan/3 = 5. For larger values of /i and M 2 , these 
bounds become clearly less stringent due to the decoupling. Notice that they are also 
quite insensitive to tan/3, since no mass insertion in Eq.(25) receives leading contributions 
from bottom Yukawa couplings. It is also worth mentioning that the bounds on the 
mass insertion (#2l)32(^rl)3i are identically to the bounds of (^ll)3i(^l)32; since they 
have the same coefficients in C* as can be seen from Eq.(26). Therefore, here we just 
present the bounds of one of them. Moreover, due to the results in Eqs.(25)-(26) and the 
strategy adopted in setting constraints, the upper bounds for the other mass insertions 
combinations, where L <-> R, turn out to be exactly the same as the corresponding ones 
in tables (2)-(4), and therefore we do not show them in our analysis. 

In analogy to the procedure used for obtaining bounds from AMs d , the imaginary 
parts will be constrained by switching off the SM CKM phase and imposing that the 
contribution of the SUSY phases to sin 2/3 does not exceed its experimental central value 
(sin2/3) cxp = 0.79. In particular we obtain 

< { Am^ ) < 28 > 

where R is defined in Eq.(26). 

In tables (5)-(7). we present our numerical results for the bounds on imaginary parts 
of mass insertions. Clearly, due to the procedure used in our analysis, these bounds turn 
out to be just proportional to the corresponding ones in tables (2)-(4), and therefore the 
same considerations about fi and tan /3 dependence hold for these bounds as well. 

Next we consider the upper bounds on the relevant mass insertions in the down-squark 
sector, mediated by gluino exchange. In Ref. [12] the maximum allowed values for the real 
and imaginary parts of the mass insertions (Sf L ) i 3 and (Sf R ) i 3 are given by taking into 
account the NLO QCD corrections. However, in that analysis the SM contributions to 
AM# d and sin 2/3 are assumed not vanishing. In order to compare our bounds on up- 
squark mass insertions with the corresponding ones in the down-squark sector, we should 
use for these last ones the same strategy adopted above. Therefore, in order to find 
conservative upper bounds on down-squark mass insertions, we will impose that the pure 
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m 


V|Im[(^J 3 i] 2 | 


VIMWsi] 2 ! 


^/\lm[(St L ) 31 (St L ) 32 }\ 


200 


1.6X10- 1 


5.4X10- 1 


2.4 xlO" 1 


400 


2.0X10- 1 


1.0 


3.0X10- 1 


600 


2.5X10- 1 


1.7 


3.8X10- 1 


800 


3.1X10" 1 


2.7 


4.6X10- 1 



Table 6: Upper bounds on mass insertions as in table (5), for M 2 = fx = 200 GeV and 
tan/3 = 5. 



m 


^|Im[(^ L ) 3 i(^)3i]| 


^/ lm[(5t L ) 31 (5 u RL U 


V / Im[(^ L ) 3 i(^)32]| 


200 


2.1X10- 1 


4.5X10- 1 


8.0 xlO- 1 


400 


3.4X10- 1 


7.2X10- 1 


1.5 


600 


5-lxlO- 1 


1.1 


2.5 


800 


7.2X10" 1 


1.5 


4.0 



Table 7: Upper bounds on mass insertions as in table (5), for M 2 = fx = 200 GeV and 
tan (3 = 5. 

gluino contribution does not exceed the experimental values on AMs d and sin 2/5, setting 
to zero the SM contribution. In these results we include the NLO QCD corrections for 
Wilson coefficients given in (22). 

The upper bounds on the real parts of relevant combinations of mass insertions {6^3)13 
(with A,B = (L,R)) from the gluino contribution to AM# d are presented in table (8). 
In table (9) we show the corresponding bounds for the imaginary parts obtained from the 
gluino contribution to CP asymmetry aj/^ Ks , again assuming zero SM contribution. The 
upper bounds on the other mass insertions in which L <-> R are not shown here, since, 
as for the chargino case, they turn to be exactly the same as the corresponding ones in 
tables (8) and (9). 

4 Light stop scenario 

In this section we will provide analytical and numerical results for the bounds on mass 
insertions, in the particular case in which one of the eigenvalues of the up-squark mass 
matrix is much lighter than the other (almost degenerates) ones. This scenario appears in 
the specific model that we will analyze in section (5.3), where the mass of the stop-right 
(rri| ) is lighter than the other diagonal terms in the up-squark mass matrix. Then the 
analytical results for the Wilson coefficients provided in section (3) will be generalized by 
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M 3 


r i 

J |Re [(S d LL ) 2 31 ] | 

• 


r i 

Vllte (^)li 1 

• 


r i 

J\Re [(5 d LL ) 31 (5 d RR ) 31 ]\ 

• 


r i 

J|Re[(^) 3 i(^ L ) 3 i]| 

• 


200 


4.6xl0" 2 


2.2xl0" 2 


8.4xl0~ 3 


1.1x10-2 


400 


l.OxlO- 1 


2.4xl0" 2 


9.6xl0" 3 


1.9x10-2 


600 


4.8X10- 1 


2.9xl0- 2 


1.2xl0- 2 


3.0x10-2 


800 


2.4X10- 1 


3.4xl0" 2 


1.4x10-2 


4.4x10-2 



Table 8: Upper bounds on real parts of combinations of mass insertions (5^ B ) 3 i, with 
(A, B) = L,R, from gluino contributions to AM Bd (assuming zero SM contribution), 
evaluated at m — 400 GeV and for some values of gluino mass M 3 (in GeV). 



M 3 


^|Im [(5 d LL ) 2 31 ]\ 


^llm [(S d RL )li]\ 


y/\1m [[S d LL ) 31 [6 d RR ) 31 ]\ 


^llm [(6 d LR ) 31 (5 d RL ) 31 ]\ 


200 


5.2x10-2 


2.5x10-2 


9.6xl0- 3 


1.2x10-2 


400 


1.2X10- 1 


2.7x10-2 


l.lxlO" 2 


2.2x10-2 


600 


5.5X10- 1 


3.3x10-2 


1.3x10-2 


3.4x10-2 


800 


2.8X10" 1 


3.9x10-2 


1.6x10-2 


5.0x10-2 



Table 9: Upper bounds on imaginary parts of combinations mass insertions (<5ab) 3 i, with 
(A,B) = L,R, from gluino contributions to sin 2(3 (assuming zero SM contribution), 
evaluated at m — 400 GeV and for some values of gluino mass M 3 (in GeV) . 

including this effect. * In our case, this modification will affect only the expression for the 
Wilson coefficient C\ [Ms) m Eq.(19), since the stop-right does not contribute to C 3 [Ms) 
at 0(A) order, as it can be seen from Eq.(20). 

By taking different the mass of the stop-right from the average squark mass, we obtain 
the following result* 

° ms) = 768^^ {l^l 2 ^ 1 ' 2 (^31 + 2A(^) 31 (^J 32 ) L 2 {x hXj ) 

- 2Y t \V il \ 2 V jl V; 2 ((5l L ) 31 (5 u RL ) 31 + X(S U LL ) 32 (5 U RL ) 31 + X(5 u LL ) 31 (8 u RL ) 32 )R 2 ( Xi , Xj ,z) 
+ YfVaV^V^^l, + 2\(5 U RL ) 31 (5 U RL ) 32 )R 2 ( Xl , Xj , z) }, (29) 

where Xi = m^+frh 2 , z = m 2 Jm 2 and the functions R 2 ( X ,y,z) and R 2 ( X ,y,z) are given 

*We do not consider here the contributions of a light right-stop to the Qi operators, since in this 
case the effect of two mass insertions (A^ fl ) 3 i can invalidate the MIA method, being no heavy squarks 
running in the loop. 

*We have used the same method introduced in Ref.[18], but our results are presented in a different 
way. 
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by 



R 2 (x,y,z) 



H 2 (x,z) 




+ Qx 2 (-1 + zf \og{x) - 6(-l + xfz 2 log(z)} 



H 2 (x,z) 



D 2 (x, z 



{(-1 + x)(x - z)(-l + z)(x + (-2 + x)z) 



+ 6x 2 (-l + ^) 3 log(a;) -6(-l + a;) 2 z(-2a; + z + log(^)} (30) 



where D 2 (x, z) — (— 1 + x) 3 (rr — z) (— 1 + z) 3 and D 2 (x, z) = (—1 + x) 2 (x — z) 2 (—1 + z) 3 
Notice that in the limit z — > 1, both the functions R 2 (x,y,z) and R 2 (x,y,z) tend to 
L 2 (x,y), recovering the result in Eq.(19). Analogously, the expressions for A 2 and A 3 
entering in Eq.(26) must be substituted by 



A 2 = Y t ^2 \Vn\ 2 VjiV* 2 R 2 (xi, Xj, z), A 3 = Y^VuV^V^R^x^z) (31) 



while Ai and A 4 remain the same. In tables (10) and (11) we show our results, analo- 
gous to the ones in tables (3)-(6), for the bounds on real and imaginary parts on mass 
insertions respectively, by taking into account a light stop-right mass. We considered two 
representative cases of m tR = 100, 200 GeV. Clearly, the light stop-right effect does not 
affect bounds on mass insertions containing LL interactions. From these results we could 
see that the effect of taking m tR < fh is sizable. In particular, on the bounds of the mass 
insertions (5^)31(^^)31, (i = 1,2) which are the most sensitive to a light stop-right. 

From the results in tables (10) and (11), it is remarkable to notice that, in the limit of 
very heavy squark masses but with fixed right stop and chargino masses, the bounds on 
(5rl)3i(5rl)3« tend to constant values. This is indeed an interesting property which shows 
a particular non-decoupling effect of supersymmetry when two light right-stop run inside 
the diagrams in Fig. (1). This feature is related to the infrared singularity of the loop 
function R 2 (x, x,z) in the limit z — > 0. In particular, we find that lim z ^ R 2 (x,x,z) = 
f(x)/x, where x = rri^/rh 2 , and f(x) is a non-singular and non-null function in x = 0. 
Then, in the limit rh >> m x the rescaling factor 1/m 2 in C* will be canceled by the 1/x 
dependence in the loop function and replaced by l/m 2 times a constant factor. 

This is a quite interesting result, since it shows that in the case of light right stop 
and charginos masses, in comparison to the other squark masses, the SUSY contribution 
(mediated by charginos) to the AB = 2 processes might not decouple and could be sizable, 
provided that the mass insertions {5^ L ) 3 i are large enough. This effect could be achieved, 
for instance, in supersymmetric models with non-universal soft breaking terms. 
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m 


™t R 


V ^ [(5l L )l] 


^/ Re[(St L ) 31 (S- RL ) 3i ] 


v /Re[(^ L ) 3 i(^ L )32] 


400 


100 


1.9X10- 1 


1.6(3.3)xl0 _1 


2.8X10" 1 


600 


100 


1.8X10- 1 


1.9(4.0)xl0" 1 


2.6X10" 1 


800 


100 


1.8X10- 1 


2.3(4.9)xl0" 1 


2.6X10- 1 


400 


200 


3.5X10- 1 


2.0(4.2)xl0~ 1 


5.2X10" 1 


600 


200 


3.3X10- 1 


2.3(5.0)xl0" 1 


4.9X10- 1 


800 


200 


3.2X10" 1 


2.8(5.9)xl0- 1 


4.8X10" 1 



Table 10: Upper bounds on real parts of mass insertions as in tables (3)-(4), for some 
values of m and fh tR (in GeV). In the fourth column the first number and the one in 
parenthesis correspond to % — 1 and % — 2 respectively. Upper bounds on mass insertions 
involving only LL interactions are the same as in tables (3)-(4). 



m 


fh tR 


>/lM0%t)§i]l 


^/|Ini[(^ L )3i(«3i]| 


V / ]Im[(^ L ) 3 i(^ L ) 3 2]| 


400 


100 


2.1X10- 1 


1.8(3.7)xl0 _1 


3.1X10- 1 


600 


100 


2.0X10" 1 


2.2(4.6)xl0- 1 


3.0X10- 1 


800 


100 


2.0X10- 1 


2.6(5.5)xl0" 1 


3.0X10- 1 


400 


200 


4.0X10- 1 


2.2(4.8)xl0- 1 


6.0X10- 1 


600 


200 


3.7X10- 1 


2.7(5.6)xl0" 1 


5.6X10- 1 


800 


200 


3.6X10" 1 


3.1(6. 7)xl0- 1 


5.4X10" 1 



Table 11: Upper bounds on imaginary parts of mass insertions as in tables (5)-(6), for 
some values of in and fh tR in (in GeV). In the fourth column the first number and the 
one in parenthesis correspond to % — 1 and % = 2 respectively. Upper bounds on mass 
insertions involving only LL interactions are the same as in tables (5)-(6). 

5 Specific supersymmetric models 

In this section we focus on three specific supersymmetric models and study the impact 
of the constraints derived in previous sections on their predictions. We discuss first 
SUSY models with minimal flavor violation, then we study the ones with Hermitian 
flavor structure, and finally we consider a SUSY model with small CP violating phases 
with universal strength of Yukawa couplings. 

5.1 SUSY models with minimal flavor violation 

In supersymmetric models with minimal flavor violation (MFV) the CKM matrix is the 
only source of flavor violation. In the framework of MSSM (with R parity conserved) we 
consider a minimal model, like in the supergravity scenario, where the soft SUSY breaking 
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terms is assumed to be universal at grand unification scale, , i.e., the soft scalar masses, 
gaugino masses and trilinear and bilinear couplings are given by 

m 2 = m 2 , M a = m 1/2 e- iaM , A a = A e~ iaA , B — B e~ iaB . (32) 

As mentioned in the introduction, only two of the above phases are independent, and can 
be chosen as 

A = arg(A*M), B = arg(£*M). (33) 

The main constraints on <pA and <pB are due to the EDM of the electron, neutron and 
mercury atom. The present experimental bound on EDMs implies that <Pa,b should be 
< 10 -2 unless the SUSY masses are unnaturally heavy [3]. 

In these scenarios, where SUSY phases <Pa,b are constrained to be very small by 
EDMs bounds, the super symmetric contributions to CP violating phenomena in K and 
B mesons do not generate any sizable deviation from the SM prediction. We have to 
mention that the universal structure for the soft breaking terms, specially the universality 
of the trilinear couplings, is a very strong assumption. Indeed, in the light of recent works 
on SUSY breaking in string theories, the soft breaking sector at GUT scale is generally 
found to be non-universal [19]. Notice that, even if we start with universal soft breaking 
terms at GUT scale, some off diagonal terms in the squark mass matrices are induced 
at electroweak (EW) scale by Yukawas interactions through the renormalization group 
equation (RGE) evolution. Therefore these off-diagonal entries are suppressed by the 
smallness of the CKM angles and/or the smallness of the Yukawa couplings. 

It is important to stress that even though one ignores the bounds from the EDMs and 
allows larger values (of order 0(1) ) for the SUSY phases 4>a,b, this class of models with 
MFV can not generate any large contribution to Ek and e'/e. Therefore, the Yukawa 
couplings remain the main source of CP violation [20]. 

Here we also found that, within MFV scenarios, the SUSY contributions to AMg d and 
dj/ipKs are negligible. In fact, due to the universality assumption of soft SUSY breaking 
terms, it turns out that the gluino and chargino contributions are quite suppressed. For 
instance, for m ~ mi/2 ~ A ~ 200 GeV and <Pa,b ~ 7t/2 (which corresponds to m 2 and 
m g at SUSY scale of order 500 GeV) we find the following values of the relevant mass 
insertions: lm(Sf 3 ) LL ~ Re(5f 3 ) LL ~ 10 -4 and lm(5f 3 ) LR ~ Re(Sf 3 ) LR ~ 10~ 6 , which are 
clearly much smaller than the corresponding bounds mentioned in the previous section. § 

Therefore, we conclude that SUSY models with MFV do not give any genuine contri- 
bution to the CP violating and flavor changing processes in K and B systems and this 
scenario can not be distinguished from the SM model one. 

5 In our analysis we have taken into account the effect of the CP violating phases in the RGE evolution. 
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5.2 SUSY models with Hermitian flavor structure 



As discussed in the introduction, a possible solution for suppressing the EDMs in SUSY 
model is to have Hermitian flavor structures [5]. In this class of models, the flavor blind 
quantities, such as the //-terms and gaugino masses, are real while the Yukawa couplings 
and A-terms are Hermitian, i.e Y^ d = Y U;d and Al d = A U)d . It has been shown that these 
models are free from the EDM constraints and the off-diagonal phases lead to significant 
contribution to the observed CP violation in the kaon system, in particular to e'/e [5]. 

Let us consider, for instance, the case of Hermitian and hierarchical quark mass ma- 
trices with three zeros [21] 



/ 



Mi = 



V 







\ 



i — u,d 



(34) 



/ 



with Ai = (m c ,m s ), B { = (m t - m u ,m b - m d ), a, = (y/m u m c , y/m d m s ), and b { = 
{^/m u m tl ^rridmiy). The phases and $ satisfy: a d — a u = n/2 and (5 d — f3 u = ir/2. 
These matrices reproduce the correct values for quark masses and CKM matrix. We also 
assume the following Hermitian A-terms: 



A d — A., — 



A 11 

A 12 e' 1 ^ 



A 12 e^ 12 
A 22 



A 13 e^ 13 "\ 
A 23 e^ 23 



\ A l3 e-^ 13 A 23 e"^ 23 A 33 J 



(35) 



Notice that, the scenario with non-degenerate A-terms is an interesting possibility for 
enhancing the SUSY contributions to Ek and e'/e [9] and it is also well motivated by 
many string inspired models. In this case, the mass insertions are given by 



m 

q 



(v«YfvS) tj v m -vY t %Ml) 



(36) 
(37) 



where q = u,d and (Y^)^ = Y^Ajj. Since the Yukawa are Hermitian matrices, they are 
diagonalized by only one unitary transformation. 

In this class of models, we find that in most of the parameter space the chargino 
gives the dominant contribution to B d — B d mixing and CP asymmetry aj/^ Ks , while 
the gluino one is sub-leading. As we emphasized above, in order to have a significant 
gluino contribution for rh ~ m g ~ 500 GeV (i.e., m ~ M±/ 2 ~ 200 at GUT scale), 
the real and imaginary parts of mass insertion {8f^) LL or (5f 3 ) LR should be of order 10 _1 
and 10 -2 respectively. However, with the above hierarchical Yukawas we find that these 
mass insertions are two orders of magnitude below the required values so that the gluino 
contributions are very small. 
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Concerning the chargino amplitude to the CP asymmetry clj/^k S i we find that the 
mass insertions {S^rl and (S^ll give the leading contribution to aj/^K s - However, 
for the representative case of m = = 200 and 4>ij — tt/2 the values of these mass 
insertions are given by 

V/Ilm^jsi] 2 ! = 6xl0" 4 , (38) 
^/|Im[(^ L ) 3 i] 2 | = 4xl0~ 3 , (39) 
v /|Im[(^ L ) 3 ife)32]| = IxlO- 4 . (40) 

These results show that, also for this class of models, SUSY contributions cannot give 
sizable effects to aj/^ Ks . As expected, with hierarchical Yukawa couplings (where the 
mixing between different generations is very small), the SUSY contributions to the B — B 
mixing and the CP asymmetry of Bd — > J/ipKs are sub-dominant and the SM should 
give the dominant contribution. 



5.3 SUSY model with universal strength of Yukawa couplings 

Supersymmetric models with small CP violating phases is a possible solution for suppress- 
ing the EDMs. In Ref. [8] it was shown that, among this class of models, the ones with 
universal strength of Yukawa couplings naturally provide very small CP violating phases. 
However due to the large mixing between different generations, it was found that the LL 
mass insertions can give sizable effects to Ek and e'/e by means of gluino and chargino 
exchanges respectively. Furthermore, it was also emphasized that in these models the 
SM contribution to the CP asymmetry aj/^ Ks might be negligible, leaving the dominant 
SUSY effect (due to the chargino exchange) to account for the experimental results. 

Here we will discuss the different contributions to Bd — Bd and aj/^ Ks in terms of 
mass insertions and compare the predictions of this model with the corresponding ones 
of Hermitian flavor structure discussed in the previous subsection. In the framework of 
the universal strength Yukawa couplings, the quark Yukawa couplings can be written as 

U,j = ^ exp and D tJ = ^ exp , (41) 

where \ u> d are overall real constants, and are pure phase matrices which are con- 
strained to be very small by the hierarchy of the quark masses [8]. The values of these 
parameters, that lead to the correct quark spectrum and mixing, can be found in Ref. [8] . 
As explained in that paper, an important feature of this model is the presence of a large 
mixing between the first and third generation. As we will show in the following, this 
property will account for large SUSY contributions to aj/^Ks- 

In the framework of universal strength Yukawa couplings Eq.(41), due to the large 
generation mixing, the EDMs impose severe constraints on the parameter space and force 
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the trilinear couplings to take particular patterns as the factorizable matrix form [8], i.e., 

( a a a ^ 

A = m b b b . (42) 

^ C C C j 

In order to satisfy the bound of the mercury EDM, the phases of the entries a, b and c 
should be of order 1CT 2 — 1CT 1 [8]. As an illustrative example, we consider m = = 
200 GeV, 4> a = 4>c = 0,0b = 0.1 and a = —1,6 = — 2, c = —3. In this case one finds 
that at low energy the average squark mass is of order 500 GeV, however one of the stop 
masses (t R ) is much lighter, mi ~ 200 GeV. The gaugino mass M 2 is of order 170 GeV 
and, from the EW breaking condition, \fj,\ turns out to be of order of 400 GeV. In this 
case, the relevant mass insertions for the gluino contribution are given by 

(5 d LL ) 31 ~ -0.001 + 0.02 i, (4l)si ^ 0.00002 + 0.0009 i. (43) 

Regarding the other mass insertions (LR and RR), they are much smaller ( < (9(10~ 6 )), 
and so we do not show them. It is clear that, with these values for the down-squark mass 
insertions, the gluino contribution to aj/^,K s is negligible (of the order of 3%). 

On the contrary, the relevant up-squark mass insertions for the chargino contribution 
are given by 

(8 U LL ) 31 ~ 0.001 + 0.05 ^, (<^ L ) 3 i ^ -0.0004 + 0.13 i, (44) 

(5 U LL ) 32 ~ -0.008 - 0.11 i, (5 U RL ) 32 ~ 0.01 - 0.28 i. (45) 

Comparing these results with the ones in tables (10) and (11), we see that for this model, 
the chargino contribution to the imaginary parts (S RL ) 3 i and (S RL ) 32 is of the same order 
of the corresponding upper bounds. Notice that these imaginary parts are of the same 
order, so that they might coherently contribute to give a sizable effect on a.j/^K s - I* 1 
particular, by using the exact 1-loop calculation, we find that the chargino contribution 
leads to sin (29 A ~ 0.75. Moreover, as a check on our computations, we have compared 
our results from MIA approximation with the corresponding ones obtained by using the 
full calculation [8]. In this case we find that, by taking into account the effect of a light 
stop, the MIA predictions are quite compatible with the results of the full computation. 



6 Conclusions 

In this paper we have studied the chargino contributions to B d — B d mixing and CP 
asymmetry aj/^ Ks in the mass insertion approximation. In our analysis we have taken into 
account the NLO QCD corrections to the effective Hamiltonian for AB = 2 transitions 
H^ =2 . We provided analytical results for the chargino contribution to H^ =2 in the 
framework of mass insertion method, and given the expressions for the B d — B d and CP 
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asymmetry clj/^k s at NLO in QCD, as a function of mass insertions in the up-squark 
sector. We have also provided model independent upper bounds on mass insertions by 
requiring that the pure chargino contribution does not exceed the experimental values 
of B — B mixing and CP asymmetry aj/^K s - Since in many SUSY models the chargino 
contribution gives the dominant effect to B — B mixing and CP asymmetry aj/^K S i our 
results are particularly useful for a ready check of the viability of these models. Moreover 
we generalized our results by including the case of a light right-stop scenario. In this 
case we found the interesting property that the bounds on mass insertions combinations 
(^l)3i(^rl)3« are n °t sensitive to the common squark mass when this is very large in 
comparison to the chargino and stop-right ones. 

Finally, we applied these results to a general class of SUSY models which are partic- 
ularly suitable to solve the SUSY CP problem, namely the SUSY models with minimal 
flavour violations, hermitian flavour structure, and small CP violating phases with uni- 
versal strength Yukawa couplings. We have shown that in SUSY models with minimal 
flavor violation and with Hermitian (and hierarchical) Yukawa couplings and A-terms, the 
SUSY contributions to the B — B mixing and the CP asymmetry cl^k s are very small and 
the SM contribution in these classes of models should give the dominant effect. On the 
contrary, in the case of SUSY scenarios with large mixing between different generations in 
the soft terms, the SUSY contributions become significant and can even be the dominant 
source for saturating the experimental value of a^x s ■ Among this class of models, we have 
investigated a SUSY model with universal strength of Yukawa couplings. In this case, we 
have found that the chargino exchange provides the leading contribution to a^x s through 
the mass insertions {b u LL )zi{b u RL )zi, i = 1,2 and (5% L ) 31 (5% L ) 32 . 
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